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ABSTRACT 

Integral field spectroscopy (IFS) represents a powerful technique for the detection and characteri- 
zation of extrasolar planets through high contrast imaging, since it allows to obtain simultaneously a 
large number of monochromatic images. These can be used to calibrate and then to reduce the impact 
of speckles, once their chromatic dependence is taken into account. The main concern in designing 
integral field spectrographs for high contrast imaging is the impact of the diffraction effects and the 
non-common path aberrations together with an efficient use of the detector pixels. We focus our 
attention on integral field spectrographs based on lcnslet-arrays, discussing the main features of these 
designs: the conditions of appropriate spatial and spectral sampling of the resulting spectrograph's 
slit functions and their related cross-talk terms when the system works at the diffraction limit. We 
present a new scheme for the integral field unit (IFU) based on a dual-lenslet device (BIGRE), that 
solves some of the problems related to the classical TIGER design when used for such applications. 
We show that BIGRE provides much lower cross-talk signals than TIGER, allowing a more efficient 
use of the detector pixels and a considerable saving of the overall cost of a lenslet-based integral field 
spectrograph. 

Subject headings: instrumentation: spectrographs — planetary system — techniques: high angular 
resolution — methods: analytical — methods: numerical 



1. INTRODUCTION 

Imaging of a significant number of extrasolar planets 
requires achieving star vs. planet contrasts of ~ 10 6 
(young giant planets), or even 10 8 - 10 10 (old giant and 
rocky planets) at a few tenths of an arc-second from a 
star, which value is ~ 10 ■ X/D in the near infrared for 
telescopes having pupil sizes of I? ~ 10 m. 

In this regime, the dominant noise contribution is due 
to the stellar background. To achieve these ambitious 
goals, high contrast imagers usually include various com- 
ponents. First, an extreme adaptive optics (XAO) sys- 
tem is used, allowing to correct aberrations up to a high 
order, and providing a high Strehl Ratio (SR). Second, 
some coronagraph is included, attenuating the coherent 
diffraction pattern of the on-axis point spread function 
(PSF). Proper combination of these two devices allows 
reduction of the stellar background down to values of 
~ 10~ a out to the AO system control radius 1 , for state- 

1 ex l/2d, d = actuator spacing projected on the telescope pupil. 



of-the-art system. This background is due to a rapidly 
changing halo of speckles generated by residual telescope 
pupil phase distortions, that have spacial frequencies 
close to those of planet images. In order to avoid false 
alarms, the detection threshold level should then be set 
at several times the root mean square (RMS) noise level. 

Even in the favourable case where the speckle in- 
tensity distribution can be assumed to be Gaussian 
(M aroi s et all l2008af ) the detection confidence limit 
should be at least 5 times the noise level. This im- 
plies that at angular separations < 10 • X/D, the lim- 
iting contrast provided by state-of-the-art extreme-AO 
and coronagraphy is ~ 2 • 10 4 for 8-10 m telescopes. In 
addition, phase aberrations originating inside the optical 
train not corrected by the extreme-AO system produce 
speckles of longer lifetime (minutes or hours) than those 
due the atmosphere. Other slowly varying (of the order 
of seconds) phase errors are due to aliasing effec ts in the 
wavefront sensor (jPoyneer &: Macintosh I l2004f ) and - 
for coronagraphic systems — to adaptive optics time-lag 
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(|Macintosh et al.ll2005f) . 

Beyond a handful of favorable case s where plan - 
ets are warm - e.g. IChauvin et al.l (|2004 I2005D : 
iNeuhaeuser et al.l (|2005l) - or with larg e separation 
from their parent star (|Kalas et al.l I2008T). or eventu- 
ally with both thes e properties (lLafreniere et alJ 120081 : 
iMarois et al.ll2008bh . additional techniques are required 
to reach the larger contrasts needed for extrasolar planets 
detection. 

Simultaneous differential imaging (SDI) is a high- 
contrast imaging differential technique by which sub- 
traction of different images of the same field acquired 
simultaneously by the same instrument allows to remove 
or reduce the noise produced by atmospheric and in- 
strumental phase aberrations. The SDI principle can 
be applied to images obtaine d with different polariza- 
tion modes (jGisler et al.ll2004 ) or se lecting two distinct 
wavelengths in a fix ed spectral range (|Lenzen et~all[2005l : 
IMarois et al.l f2005^) . or better exploiting th e entire spec- 
tral r ange by integral field spectroscopy (jBerton et al.l 
2006). In this paper we will focus on SDI based upon 
this latter strategy only. 

Essen t ially, SDI is a ca li bration technique 
(ISmith I Il987t iRacine et all 119991: IMarois et all 120001: 
Sparks fe Fordl 12002 1 : iBiller et all l2004t iBerton et all 
20061: iRen fe Wang 1 120061 ; iThatte et al.M2007l) : images 
are acquired simultaneously in bands at close wave- 
lengths where the planetary (but not the stellar) flux 
differ appreciably. Subtracting each other these images 
should allow to remove or at least reduce the speckle 
noise, since this is assumed to be similar in the various 
images after a suitable chromatic re-scaling, while the 
planet signal is left nearly untouched. 

There at least two ways to exploit this calibration 
technique. In the more traditional approach, specific 
characteristics of the (expected) planetary spectra are 



(ISudarskv et al.ll200C ; iBaraffe et al.ll2002: 
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2004) and ob- 



servations (e.g. of Brown Dwarfs and gaseous planets 
in the Solar System), the spectra of giant planets are 
dominated by several absorbtion bands (mainly due to 
methane and water vapor) at both visible and near in- 
frared wavelengths. In such a case, SDI may work by 
subtracting images where the planet signal is absent from 
those where it is present, while the background is nearly 
the same, because the spectrum of the parent star is 
nearly featureless, see Figure [T] The main advantage of 
this technique, is the minimum assumptions required on 
the chromatic behavior of speckles; however, this tech- 
nique allows only a limited reduction of noise. Alter- 
natively, we might t ry to model the varia tion of speck- 
les with wavelength (Sparks fc Fo rd 2002). In principle 
this allows to remove completely speckle noise without 
making any assumption about the planetary spectrum, 
hence allowing to r etrieve the real planetary spectrum 
(jThatte et al.|[2O07h . 

Independently to the adopted SDI recipe, integral field 
spectrograph designs tuned for diffraction-limited high- 
contrast imaging should take into account several effects 
jeopardizing the interpolation procedures requested be- 
fore simultaneous spectral subtractions, which in turn 
severely limit the accuracy of this calibration technique. 

In this paper we present a discussion of these effects 



0-8 1.0 1.2 1.4 1.6 1.8 

Wavelength (micron) 

Fig. 1. — Near infrared low resolution spectra (two-pixel resolving 
power R = 50) of the Brown Dwarf GL 570 D, of the planet Saturn 
and of a G2V star (by courtesy of Dr. L. Testi and Dr. F. Ghinassi) 
obtained at TNG+NICS with its Amici Prism. The spectra are 
normalized to their flux at 1.30 [ira. 



and derive the basic equations that should be considered 
when designing lenslet-based diffraction-limited integral 
field spectrographs. Then, we describe a new concept for 
the lcnslct-array shaping the IFU of such instruments 
(i.e. BIGRE) allowing to improve significantly over the 
main limitations of the more traditional designs based on 
the TIGER concept. 

The structure of the paper is as follows. In § [2] we re- 
call the basics of a post-coronagraphic speckle field. In 
§ we summarize the principle of SDI. In § |4] we dis- 
cuss the basics of spectroscopic SDI (hereafter S-SDI), 
defining the conditions allowing to avoid aliasing errors 
when sampling both the entrance speckle field and the 
final exit slit functions. In §[5] we present various options 
for IFS-instruments suited for S-SDI. In §[6] we define the 
cross-talk terms in the case of diffraction-limited lenslet- 
based IFS. In §[7] we derive the rules governing the image 
propagation at the diffraction limit through the TIGER 
concept, and in § [8] the ones proper to the new BIGRE 
concept. Specifically, we explain here how to conceive a 
BIGRE-orientcd IFS instrument adopting standard diop- 
tric devices. In §[9] we present two design setups (based 
on BIGRE and TIGER respectively) for SPHERE 2 , indi- 
cating the solution adopted for its future IFS. In SlTOlwe 
compare the TIGER and the BIGRE concepts in terms 
of coherent and incoherent signals suppression, consider- 
ing several cases for the single lens shape and the IFU 
lattice configuration. Finally, our conclusions are drawn 
in §HI1 

2. POST-CORONAGRAPHIC SPECKLE FIELD MODELING 

An approp riate underst a nding of chromatic inten- 
sity (e.g. I Racine et all (|1999f) ) and spatial (e.g. 
iSparks fc Fordl (|2002f ^ scaling of a speckle field is ba- 
sic to any application of the SDI calibration technique. 
For this reason a short description of these physical con- 

2 SPHERE is an instrument designed and built by a consor- 
tium of LAOG, MPIA, LAM, LESIA, LUAN, INAF, Observatoire 
de Geneve, ETH, NOVA, ONERA and ASTRON in collaboration 
with and under from ESO. Its science objective is the direct detec- 
tion and character ization of giant ext rasolar planets in the visible 
and near- infrared HBeuzit et al.ll2008h . 
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cepts is fundamental to introduce the reader to the topics 
tr eated in the rest of the paper. Inspired by the approach 
of lPerrin et aHd2003h . we will use the Fraunhofer approx- 
imation to describe the impact of small residual phase 
variations of the electric field (e) imaged on a fixed post- 
coronagraphic entrance pupil plane 3 , i.e. the working 
case of high-contrast imaging instruments like SPHERE. 

While this approach allows a simple mathematical 
treatment and physical understanding, it ignores more 
complex effects due to am plitude errors and Fr esnel prop- 
agation, as pointed out bv lMarois et a l. ( 2006). It is out- 
side the scope of this paper to discuss such effects, which 
can be minimized by careful instrument design, but it 
is likely that they will set the ultimate limit of planet 
imaging. 
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Fig. 2. — Example of a post-coronagraphic speckle field (inte- 
gration time = 0.5 msec , wavelength = 1 /im), simulating the 
extreme- AO system and the 4 quadrant phase mask coronagraph 
of SPHERE (by courtesy of the SPHERE team). 



The most general expression of the monochromatic 
electric field once projected on the coronagraphic en- 
trance pupil plane is: 

e = p ■ exp[i ■ (/>], (1) 

where (p) is the coronagraphic pupil transmission func- 
tion, and ((f>) is the phase of the electric field evaluated 
over this coronagraphic pupil plane. Assuming a perfect 
optical propagation from the telescope to this plane - i.e. 
no differential chromatical aberrations in the beam - the 
chromatism of the phase can be written explicitly as a 
function of the wavelength (A) and the wavefront error 
(w) as follows: 



2tt 

T 



(2) 



Assuming as real the expectation value of the wavefront 
error given by an extreme- AO system in the near infrared 
(i.e. w < 10 -2 fim at A ~ 1 /xm), equation ([!} can be 
approximated as follows: 



V 0- + i 



(3) 



3 Hereafter Fourier pairs are defined with the same letter written 
in small and capital case respectively. 



At this point, the action of an un-specified corona- 
graph can be formalized directly on the coronagraphic 
exit pupil plane. The goal of the coronagraph is to can- 
cel as much as possiblee the amplitude of the electric field 
along the optical axis on this plane. Exploiting ([3|), the 
resulting on-axis electric field (e c ) for a perfect corona- 
graph 4 is then: 

e c = e - p = i ■ p ■ <f>, (4) 
or, by equation @, is equal to: 

e c = l- — -p-w. (5) 

Defining finally (E C ,P,W) as the Fourier transforms 
(FT) of (e c ,p,w), equation §5§ allows to express the 
monochromatic post-coronagraphic speckle field (S) as: 

S(\)^\E C (\)\ 2 = 0^) .|P®wf . (6) 

Equation © shows that the intensity of a speckle field 
scales proportionally to A~ 2 , while its chromatic wave- 
length scaling comes from the fact that the variable in- 
volved in the wavefront w is the spatial frequency (v) 
and not the position (x) in the image plane, i.e.: 



w{v) = FT[W{x)}. 



(7) 



This indicates that spatial frequency translates into posi- 
tion according to wavelength, e.g. by applying the stan- 
dard grating equation as follows: 



m ■ A = g • sin(8), 



(8) 



where (m) is the diffraction order, (9) the diffraction an- 
gle and (g) is the grating constant corresponding to the 
spatial frequency (v), or: 



giy) 



the position on the image plane returns: 
x = f ■ sin(0), 



(9) 



(10) 



/ being the focal length of the post-coronagraphic re- 
imaging optics. Using equations © and © this position 
can be written finally as: 



x = f- 



A- 



(11) 



Equation (jlip indicates that the position of a speckle cor- 
responding to a given fixed spatial frequency due to the 
post-corona graphic wavefront err or scales linearly with 
wavelength (|Sparks fc Fordl [20021 ). More in detail, this 
means that for every fixed position in the image plane, 
speckles corresponding to distinct spatial frequencies get 
distinct wavelengths (Figure [3]). We call this feature 
speckle chromatism. 

4 A perfect coronagraph removes actually the coherent part of 
the electric field amp l itude due to the on-axis optical beam only, see 
e.g. ICavarroc et al. (2006); here we consider the total amplitude 
in order to simplify the related formalism. 
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Fig. 3. — Cuts through the data cube obtained by the chromat- 
ical dispersion of a post-coronagraphic speckle pattern resulting 
from end-to-end simulations of the IFS inside SPHERE (intensity 
increases from red to blue colors). Position on sky (9) is on the 
horizontal axis, while the spectral range (X m i„ — Xmax) is on the 
vertical axis. The red and blue lines indicate two spectra taken at 
different radial distances to the optical axis. Moving along these 
spectra, a variable pseudo-periodic modulation due to the speckle 
chromatism is clearly visible (by courtesy of the SPHERE team). 

3. THE SDI CALIBRATION TECHNIQUE FRAMEWORK 

In the approach considered in this paper, the funda- 
mental SDI step is the simultaneous acquisition of im- 
ages at adjacent wavelengths in a spectral range where 
the planetary and stellar spectra differ appreciably. From 
ground-based observations, the wavelength bands Y, J, 
H, and K are well suited for extrasolar g iant planets 
(iBeuzit et ail 120081: Macintosh et all 120081), and rocky 
planets (jVerinaud et al.ll2008l ). 

Let S(A, 9) be the monochromatic spectral signal corre- 
sponding to a fixed angular position (9) on sky expressed 
as the sum of the spectral signal of the star, Si (A, 9), and 
the spectral signal of a candidate low-mass companion 
(e.g an extrasolar planet) which lies specifically in this 
angular position, Pl(X,9). Fixing a pair of wavelengths 
(Ai, A2) inside the window above, the following relations 
hold: 

S(Ai,0) = Si(Ai,0) + PZ(Ai,0) (12) 
S(X 2 ,9) = St(X 2 ,9)+Pl(X 2 ,9). (13) 

The basic SDI assumption is that after suitable flux 
normalization and chromatic re-scaling, the following re- 
lations hold for the boundary wavelengths of the range 
above: 

St(X u 9) = St{X 2 ,9) (14) 

pi(x 2 ,e)=o. (15) 

Then the difference between S(Ai, 9) and S(A2, 9) should 
return — in principle — the spectral signal Pl(Xx,9) 
only, i.e. the one appropriate to the low-mass (or ex- 
trasolar planet) candidate. However, while working with 
narrow-band filters several precautions are required: 



• an image taken with one filter has to be spatially 
re-scaled before confronting it with a second im- 
age taken with a different filter due to the speckle 
scaling described in § [2] 

• any filter separating two adjacent spectral bands 
should have similar spectral transmission profiles; 

• the difference (§\ij) between the central wave- 
lengths (Aj,Aj) of two adjacent filters should be 
as small as possible. 

The last item is the most critical due to the fact that 
chromatism of the speckle field always induces a cer- 
ta in amount of phase errors. Adopting the formalism 
of iMarois et al.l (|2000l ). the residual wavefront distortion 
can be described through the Fourier transform of the 
post-coronagraphic wavefront error (W), or by its rela- 
tive chromatic phase-error (<J>). 

Ado pting the standa rd approximation for the Strehl 
ratio (|Marechal I Il947f) . it is possible to transfer this 
RMS wavefront-error (c$) on a relative flux variation 
on the detector plane. In detail, defining <t$(A) as 
the RMS chromatic wavefront-error- TMarois et al.l (2000) 
found the following relation for the flux residual between 
images taken with two narrow-band filters 

^i = 2 .al(X i ). S -^. (16) 

Equation (fT6|) indicates that with the so called single 
difference method the final error is proportional to: 

• the variance of the wavefront error: erj, 

• the relative wavelength separation between the 
narrow-band filters: 6Xij/Xi. 

The need for a calibration technique more efficient than 
SDI but still based on the simultaneous difference of chro- 
matic image s of the same t arget field was addressed the- 
oretically bv lMarois et a l. (2000), which showed that the 
speckle noise reduction could be much more efficient if 
observations at three wavelengths were available using 
their double difference method, and tasted experimen- 
tally with the discovery the first planet obtained by using 
this calibration technique (lLagrange et al.ll2008f) . whose 
infrared contrast is > 10 2 . 

Starting from there, it is reasonable to assume that a 
larger number of images at different wavelengths, taken 
with a regular spectral-step, can result in even better 
reduction of speckle noise with a true S-SDI calibration 
technique. The gain could be even larger if observations 
at several wavelengths would allow an accurate deriva- 
tion of the chromati c wave length scaling, as proposed 
e.g. bv iThatte et all (|2007l ). This thought suggests the 
use of integral field spectroscopy for collecting data si- 
multaneously at a large number of wavelengths given by 
the total spectral l ength and the spec tral resolution of a 
suitable disperser (|Berton et al.ll2006l ). 

Note that such an approach is convenient even in the 
more conservative approach where modeling of the spec- 
tral dependence fail, simply because a larger number of 
wavelength pairs can be constructed. 



BIGRE 



5 



4. THE S-SDI CALIBRATION TECHNIQUE FRAMEWORK 

Exploiting an IFU as field stop array over an opti- 
cal plane conjugated with the focal plane of the tele- 
scope itself allows an appropriate sampling of the post- 
coronagraphic speckle field defined by equation © . The 
fact that this optical signal gets a finite cut-off spatial 
frequency proportional to D/X m i n , where D is the post- 
coronagraphic pupil size and X m in is the spectrograph's 
cut-on wavelength, means that a correct spatial sampling 
on this plane should be imposed searching for suitable 
sizes for the separation between adjacent spaxels 5 , which 
in turn compose the adopted IFU. This sampling condi- 
tion is detailed in § 14.11 

The request of a sampling criterion based upon the 
Shannon theorem is mandatory not only at the level of 
the IFU spaxels but also at the level of the detector pix- 
els. In this case the Shannon sampling condition allows 
to interpolate correctly, both spatially and spectrally, the 
exit slit functions, which in turn are the final output of 
an integral field spectrograph. These two sampling con- 
ditions are detailed in § 14.21 and § 14.31 respectively. 



4.1. Spatial sampling of the entrance speckle field 

Let Fi n be the focal ratio by which the post- 
coronagraphic speckle field is projected on the IFU plane. 
Theory of image formation (e.g. Goodman 1996) implies 
then that the cut-off spatial frequency appropriate to S 



Fi n ■ X 



in "mm 



(17) 



The spaxel size (D spaxe i) defines the Nyquist spatial 
frequency on this plane: 



VNy 



1 



2 • D 



spaxel 



(18) 



Thus, the Shannon sampling theorem applied to the 
IFU plane returns: 



VNy > VC- 



(19) 



4.2. Spatial sampling of the spectrograph's exit slits 

The condition avoids aliasing effects when interpolat- 
ing the array of exit slits over the whole range of wave- 
lengths considered by the spectrograph, and it may be 
written through the following formalism. 

The detector pixel size {d P i xe i ) defines the Nyquist spa- 
tial frequency on this plane: 



(J-Ny 



1 



2-d 



pixel 



(20) 



Once the final spectrograph's exit slits are imaged on the 
detector pixels through a fixed output focal ratio (F out ) 
and an optical magnificat i on [m jjrg), theory of image 
formation, e.g. iGoodmaril (|1996f ). implies that their spa- 
tial cut-off frequency is: 

5 Spaxel indicates a spatial pixel appropriate to the IFU sub- 
system inside an IFS-instrument. IFU in turn is the matrix of 
spaxels which should be placed on the re-imaged telescope focal 
plane, working as an optical field-stop array. 



1 



A r , 



miFs ■ F . 



(21) 



where \ m in indicates the shortest wavelength imaged by 
the spectrograph. We define the Super-sampling condi- 
tion as: 



(22) 



4.3. Spectral sampling of the speckle field over the entire 
field of view 

When working with a speckle pattern data-cube, chro- 
matic re-sampling is needed to obtain both monochro- 
matic images, as indic ated by iMarqis et al.1 (|2000f ). or 
spectra, as in dicated bvlThatte et all (|2007h . 

To this aim lSparks fe FordT i 2002f ) suggested to adopt 
a suitable pixel-dependent re-sampling of the speckle 
field which varies according to wavelengths, while 
iRen fc Wang I (|2006h developed a subtraction algorithm 
based upon analytical modelings of the spectral content 
of a speckle field. Anyhow, before any re-sampling recipe, 
it is important to find out the exact condition allowing 
to avoid aliasing errors due to the speckle chromatism 
effect. 

Since the speckle pattern scales proportionally to wave- 
length (§ O, a feature located at an angular distance 9 
from the central star at wavelength A moves spectrally 
at a rate of dX/d9 — X/6. Spatial speckles of width 
S9s = X/D therefore translates into spectral speckles of 
width: 



6X S = 



A 2 



9-D' 



(23) 



i.e., the spectral extension of speckles is inversely pro- 
portional to the distance from the field center. Nyquist 
sampling of spectral speckles requires spectral sampling 
(SXp) corresponding to half the speckle width, so far a 
two-pixel resolving power (R = X/2SXp), Nyquist sam- 
pling implies the following condition: 



R > 



X 
SXl 



D 
A"' 



(24) 



This condition will be fulfilled within a field angle 9n v , 
referred to as the Nyquist radius, given by: 



9 Ny =R - "p- 



(25) 



We note that it is possible to ensure Nyquist sampling 
in a system which does not fulfil the Super-sampling con- 
dition written in equation (f2"2")h as long as its field of 
view does not exceed the Nyquist radius and as long as 
the source itself does not contain spectral features which 
violate the Shannon theorem. For example, an instru- 
ment operating on an 8 meter telescope at 1.6 /im with 
a full field of view of 5 arc-seconds, would require a two- 
pixel resolving power of at least 60. For systems where 
larger field of view or lower resolving power is required, 
the Super-sampling condition must be fulfilled. In these 
systems, the zone lying within the Nyquist radius fulfils 
both equation (|2"2"|) and and equation (j2~4"|) . We refer to 
this double fulfilment as Hyper-sampling. 
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For an integral field spectrograph covering a spectral 
range fixed between a cut-on (X m in) and a cut-off wave- 
length (A max ), where A c represents the central one, the 
Hyper-sampling condition will be valid over the whole 
spectral range within the radius: 

0n v = (26) 

5. OPTIONS FOR THE IFS CONCEPT REALIZING S-SDI 

IFS needs a very large number of pixels at the level of 
the final image plane where the matrix of spectra is ac- 
quired by the detector. This issue is particularly impor- 
tant when spectral and spatial information are recorded 
simultaneously in the detector plane, such as for IFS 
based on the image sheer or the TIGER 6 concepts. 

The image sheer option is more efficient in terms of 
detector pixels usage, since no separation between spec- 
tra from adjacent pixels is required in one space dimen- 
sion. Assuming a square detector, the number of detec- 
tor pixels (iVj et ) required for a given number of spax- 
els {N^ paxel ) and number of spectral samples (N spec ), is 
given by the following relation: 

Ndet — ^spaxel ' N spec . (27) 

In this concept a bi-dimensional field of view is divided 
by mirrors into strips, and then re-formatted on a mono- 
dimensional pseudo long-slit (see Figured]). Monochro- 
matic exit slits will be then obtained downstream, by 
using a standard collimator, disperser and camera op- 
tical system. A potential problem of the image sheer 
design concerns the non common path aberrations in ad- 
jacent spaxels of the field of view that fall on different 
slices. However this concept has been proved able to 
obtain (moderately) high-contrast images from ground 
even without coronagraphic devices and with m oderate 
Strehl ratios (~ 0.3 - 0.5) (|Thatte et al.ll2007ft . A fur- 
ther examination of an image sheer instrument dedicated 
to high-contrast diffraction-limited imaging spectroscopy 
is on progress within the feasi bility study for the future 
E-ELT Planet Finder facility (jKasper et al.ll2008[ ). 

On the other hand, non-common path aberrations are 
expected to be very small i n the case of the TIGER-type 
concept (|Bacon et al.lll995[) , which uses an IFU based on 
a matrix of lenses with fixed lens pitch. In this case spec- 
tra given by individual spaxels should be separated on 
the detector. For a separation of N sep between spectral 
samples, the required number of detector pixels becomes: 

N det — N spaxer N spec ' N sep . (28) 

The lenslet-based concept then requires a large number 
of detector pixels. However, the format of image sheer 
IFS-data on the detector is suited for spectra with many 
spectral elements, i.e. > 10 2 , and relatively small num- 
ber of spaxels, i.e. < 10 4 . This are not typical values 
for instruments dedicated to planet search that gener- 
ally requires short spectra (~ 20 — 30 spectral elements) 
for a large number of spaxels (~ 10 5 ). In order to ade- 
quately exploit the detector, the number of slices should 

6 "TIGER" is a French acronym standing for "Traitement 
Integr al des Galaxies par l'Etude de leurs Rays" , has Bac on et al.l 
(1995) named their lenslet-based IFS. 




Fig. 4. — The principle of the image slicer IFU, as figured in Pri- 
eto & Vives 2006 (by courtesy of Dr. Eric Prieto and Dr. Sebastien 
Vives) . 



Telescope Spectrograph Spectrograph 

focal surface input output spectrum 




Detector 

Fig. 5.— The principle the TIGER IFU as figured by Lee et 
al. 2001 (by courtesy of Dr. Jeremy Allington-Smith and his co- 
authors). 

be roughly given by the ratio between the spaxels and 
the length of the spectra. This value is ~ 10 3 for an in- 
tegral field spectrograph tuned to planet finding, which 
would result in an extremely long pseudo-slit. The for- 
mat of the image slicer IFU then exacerbates the prob- 
lems related to non common paths: indeed photons from 
adjacent spaxels may have very different paths through 
the instrument. It is then difficult to maintain small 
the phase errors, possibly compromising most demand- 
ing high-contrast imaging. 

Given the difficulties inherent to the image slicer solu- 
tion, we carefully examined the properties of the lenslet- 
based design, trying to minimize the separation between 
spectral samples. To this aim, we developed t he new op- 
tical concept proposed bv lDohlen et al.l (|2006l ): BIGRE 7 . 
The properties of this design are discussed and compared 
to the TIGER ones, starting from § [5] 

6. INCOHERENT AND COHERENT CROSS-TALKS OF A 
LENSLET-BASED IFU 

Adopting the formalism of iGoodmanl (|1996f ). any 
spaxel of an IFU is a sum of linear optical systems. In 
the specific case of a lenslet-based IFU these systems 
are the single lenses. The coherent and incoherent part 
of the electric field incoming onto these optical linear 

7 "Bigre" was the first word uttered by G. Courtes - the inventor 
of the TIGER concept - while the Authors explained him all the 
problems of diffraction-limited IFS and their possible resolution 
using this new optical concept. "Bigre" is a French exclamation 
with a meaning similar to the British: "Bligh-me" or the Italian: 
" Accidenti" . 
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systems are transmitted in a different way through two 
adjacent spaxels. Specifically, when the illumination is 
coherent, the linear responses of adjacent spaxels vary 
in unison, and therefore their signals, once transmitted 
and re-imaged on the spectrograph's slits plane, must be 
added in complex amplitude. Contrarily, when the illu- 
mination is incoherent, the linear responses of two ad- 
jacent spaxels are statistically independent. This means 
that their signals, once transmitted and re-imaged on the 
spectrograph's slit plane, must be added in intensity. 

Hence, once dispersed and re-imaged by the spectro- 
graph's optics 8 , monochromatic slits corresponding to 
adjacent spaxels will suffer from a certain amount of 
interference. We call this quantity coherent cross-talk. 
Furthermore, monochromatic slits will be affected by a 
spurious amount of signal due to its adjacent spectra. 
We call this quantity incoherent crosstalk. With refer- 
ence to Figure |6l coherent cross-talk is the interference 
signal between monochromatic spectrograph's entrance 
slits which correspond to adjacent lenses, i.e. separated 
by a distance equal to the IFU lens pitch 9 . While, in- 
coherent cross-talk is the spurious signal registered over 
a fixed monochromatic spectrograph's exit slit and due 
to its closest spectra, even if due to photons of different 
wavelength. 




Fig. 6. — Sketch of the final spectra (black rectangles) superim- 
posed on an array of 7 (red) hexagonal spaxels, down of a lenslet- 
based IFU. Dl indicates the IFU lens pitch, while the dispersion 
axis is indicated through a black arrow labeled with the symbol A. 



Incoherent and coherent cross-talks represent a major 
issue identifying the best solution for the spaxels shape 
(circular, square, etc.), the lenslet lattice configuration 
(hexagonal, square, etc.), and for the geometric allo- 
cation of the spectra at the level of the detector plane. 
In fact, incoherent and coherent cross-talks are spurious 
signals — not removed by the application of Super- and 
Hyper-sampling criteria — which still affect the final ar- 
ray of spectra, thus damaging the final three-dimensional 

8 The dispersion axis can be defined orienting the spectrograph's 
disperser with respect to a reference frame fixed on to the IFU. 

9 The pitch of an array of spaxels is defined as the center-to- 
center distance among adjacent ones. For a filling factor close to 
unity this quantity equals the size of the single spaxel. 



data cube. The selection of the kind of field unit to be 
mounted at the entrance of a lenslet-based integral field 
spectrograph should then depend on the estimate of the 
level of incoherent and coherent signals over the individ- 
ual exit slits of such a spectrograph. Additional consid- 
erations should enter in this choice, e.g. the fact that the 
relevance of the cross-talk terms depends on the wave- 
front errors after the coronagraph or that minimization 
of the cross-talk might result in a system design which is 
potentially less efficient when observations are limited by 
photon noise. In general, cross-talk should be specified 
so that its contribution to the contrast error budget is 
less than the flat field errors and all remaining spurious 
effects affecting the the post-coronagraphic speckle field. 

6.1. Coherent cross-talk: the formalism 

Basically, coherent cross-talk is the interference of a 
beam passing through a number of apertures (individual 
lenslets) and measured on a screen (the spectrograph's 
entrance slits plane) conjugated to the detector plane. 

Let us assume a flat wavefront impinging onto the IFU 
lenses. Let now be E\ the complex electric field of the 
coherent signal transmitted by spaxel 1 on the spectro- 
graph's entrance slits plane. Let be dE 2 the stray part 
of the complex electric field of the coherent signal trans- 
mitted by spaxel 2 (spaxel 2 being adjacent to spaxel 1) 
and evaluated in the position of the slit corresponding to 
spaxel 1. Ei and 0IE2 are complex quantities that dif- 
fer according to the phase difference, which is due to 
different optical paths through the different apertures 
(lenslets). The effective coherent intensity measured on 
the spectrograph's entrance slits plane and correspond- 
ing to the position of spaxel 1 will then be: 

if = \E X + dE 2 \ 2 . (29) 

In the worst case, the phase difference of waves passing 
through adjacent lenses is irk (k € Z). In this case and 
neglecting the term \dE 2 \ in the binomial expression of 
equation (|29p . the effective coherent intensity proper to 
spaxel 1 becomes: 

Jf = h + 2 • 1^1 • \dE 2 \ = h ■ (1 + CCT). (30) 
CCT is defined as the coherent cross-talk coefficient: 

where the stray coherent intensity proper to spaxel 2 
evaluated in the position of spaxel 1 is defined as: 

dI=\dE 2 \ 2 , (32) 
and the own coherent intensity of spaxel 1 is defined as: 

/i = |£i| 2 (33) 

CCT represents the maximum extra-amount of coher- 
ent signal on the slit function corresponding to a fixed 
lenslet aperture, and its estimate can be given by mea- 
suring the square root of the coherent intensity proper 
to the slit function corresponding to the adjacent aper- 
ture. However, the total amount of coherent cross-talk is 
obtained only by adding the contribution due to all the 
apertures in the lenslet-array. 
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6.2. Incoherent cross-talk: the formalism 

The amount of spurious incoherent light can be eval- 
uated directly on the detector plane, where a single exit 
slit appears as a spectrum. As indicated in Figure [5J any 
final spectrum is surrounded by several adjacent spectra. 

Let be I\ (A) the intensity proper to a fixed monochro- 
matic exit slit; due to the presence of an adjacent exit 
slit its effective incoherent intensity will be: 

l{ NC {\) = h{\) + dh{X), (34) 

where dl2{X) is the stray incoherent monochromatic in- 
tensity of a given adjacent exit slit, evaluated at a dis- 
tance equal to the separation to the fixed one: 

dI 2 (X) = ICT(X) ■ h(X) (35) 

where ICT(X) is defined as the monochromatic term of 
the incoherent cross-talk coefficient (ICT). 

The incoherent cross-talk coefficient corresponding to 
the spectrograph's wavelengths range (A m j n — X max ) is 
then defined as: 

Thus — differently to the coherent case — the incoher- 
ent cross-talk must be considered on the detector plane, 
searching for spectral alignments for which the distance 
among adjacent spectra is minimized. Once this spectral 
alignment is found, an estimate of ICT can be given by 
measuring the incoherent intensity of a single monochro- 
matic exit slit at the distance equal to the transversal 
separation among adjacent spectra. However, the to- 
tal amount of incoherent cross-talk is obtained only by 
adding the contribution of all the spectra imaged onto 
the detector plane. 

7. DIFFRACTION-LIMITED INTEGRAL FIELD 
SPECTROSCOPY WITH THE TIGER CONCEPT 

In classical TIGER design optimized for seeing lim- 
ited conditions the spaxels (or microlenses) composing 
the IFU are much bigger than the Airy disk, providing 
therefore resolved images of the telescope entrance pupil, 
which in turn represent the entran ce slits of thi s kind 
of int egral field instrument, see e.g. iBacon et al.l (|1995l . 
[200l . 

Differently, in the case of high-contrast imaging the 
microlenses sample the telescope image according to the 
Shannon theorem. Each microlens acts like a diaphragm 
isolating a portion of the incoming electric field and con- 
centrates it into a micropupil image in the focal plane 
of the microlens, acting as the entrance slit function of 
the spectrograph. The micropupil image is the convolu- 
tion between the geometrical pupil image and the PSF 
of the microlens. As seen below, Nyquist sampling of the 
focal plane implies that the telescope entrance pupil is 
unresolved by the microlens. 

For a circular lens of diameter D spaxe i the transmis- 
sion function is H(u / D spa xei) , where u is the image co- 
ordinate normalized to the lens diameter and H(x) is a 
top-hat function with unitary transmission within the 
unitary diameter and zero outside of this diameter. Ac- 
cording to condition (Til?)) , the size of the single microlens 
should be: 



D spaxel < 2 A """ ) . (37) 

Following Born & Wolf] (|I965| ), the monochromatic full- 
width-half-maximum (FWHM) of the PSF proper to a 
circular microlens with focal length f out is: 

FWHM = f out ■ 1 .02 • — (38) 

D spaxel 

while the geometrical diameter of the micropupil is: 

D MPG = (39) 

"in 

Combining equations (|37| , (|38|) and (|39|) we obtain then: 
FWHM = 2.04 • ( Dmpg ' X \ , (40) 

This size is therefore at least twice as wide as the geo- 
metrical pupil, and so the convolution product is approx- 
imately equal to the microlens PSF. 

Thus, we can stay that the field distribution onto the 
spectrograph's slit plane approximates the one proper to 
an unresolved micropupil, which is described by the Jinc 
function 10 corresponding to the microlens aperture: 

epu P ii(s) ~ Jinc(s), (41) 

s being defined as the pupil the co-ordinate normalized 
to A • F^, where F out is: 

F ut = (42) 

-LJspaxel 

Finally, the slit function will be the square modulus of 
this signal: 

SF(s) = \e pupU {s)\ 2 . (43) 




Fig. 7. — Normalized TIGER entrance slit function: the example 
shows the case of the IFU optimized for SPHERE in the working 
wavelengths range: 0.95 — 1.35 /im; colors indicate seven distinct 
wavelengths. 



10 We define Jinc function as the Fourier transform of circular 
aperture: Jinc(x) = (2 ■ Ji(tt ■ x)) / (n ■ x), where J\ indicates the 
Bessel-J function of order one. 
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7.1. Sampling analysis applied to the TIGER concept 



As indicated by equation (j4"5|) , the single spectro- 
graph's slit is an un-bound signal whose size varies lin- 
early with wavelength. The final pixel size defines the 
spatial Nyquist frequency on the spectrograph image 
plane according equation (|20|) . Due to its un-bound na- 
ture, the spatial cut-off frequency of the spectrograph's 
exit slit gets the finite value fixed by equation (pH]) . Then, 
following condition (|22|) , Super-sampling imposes a lower 
limit to the output focal ratio by which the single mi- 
crolens generates its corresponding micropupil: 



F n 



> 



2-d 



pixel 



An 



miFS J 



(44) 



Output focal ratios lower than the one fixed by equa- 
tion (|44|) introduce aliasing errors in the sampled spec- 
trum, unless the field is smaller than the Nyquist radius. 
According to condition (j2"6")l , this latter depends on the 
post-coronagraphic pupil size, the spectrograph's work- 
ing wavelengths range and its spectral resolution. Hence, 
the true Hyper-sampling is obtained when this radius 
matches with the maximum image field radius, which 
in turn is related to the spectrograph's resolving power. 
Then, for a fixed resolving power, Hyper-sampling is then 
a matter of allocation of the array of final spectra onto 
the detector pixels, which in turn depends on the ac- 
cepted cross-talk levels. 



TABLE 1 

Independent parameters of a TIGER-oriented IFS 

Post-coronagraphic pupil size = D 

IFS cut-on wavelength = A m i n 

Size of the single TIGER microlens = D apaxe [ 

Focal length of the TIGER microlens = f ou t 

IFS detector pixel size = d pixe i 

IFS optical magnification = m/pj 

IFS disperser (2-pixel) resolving power = R 



8. DIFFRACTION-LIMITED INTEGRAL FIELD 
SPECTROSCOPY WITH THE BIGRE CONCEPT 

Cross-talk in diffraction-limited TIGER-oriented IFU 
is generally quite large because the output slit functions, 
taking the form of an Airy pattern, decrease slowly with 
distance from the center. Suitable apodization of the mi- 
crolenses might in principle be used to reduce the cross- 
talk terms, but the feasibility of such a scheme remains 
to be demonstrated. We consider instead an alternative 
lenslet-based optical scheme that we call BIGRE, which 
we consider to be much more practical. 

As in the TIGER case, the BIGRE spaxel consist of a 
microlens which acts essentially as a diaphragm isolating 
a portion of the incident electric field. This lens, of focal 
length fx, focalizes the field into an unresolved microp- 
upil with a field distribution described by equation (|41l) . 
Differently to the TIGER case, we place a second mi- 
crolens at a distance equal to its focal length /2, behind 
the micropupil. This lens collects field and reproduces 
an image of the first lens, behind the micropupil. When 
/2 < /i, the final image is reduced, resulting in the same 
flux-concentrating effect as in the original TIGER con- 
cept, but without the field-pupil inversion. We define K 
factor as the spaxel de-magnification factor: 



J2 



(45) 



Ideally, for infinitely wide optics throughout the follow- 
ing spectrograph, the slit function is a perfectly bound 
top-hat function, so no cross-talk would be present be- 
tween spaxels. Of course, this is not physical, and the 
following finite sized optics modifies the slit function as 
we will see in the following. 

It may also be argued that a perfect top-hat function 
is not the ideal slit function from a sampling point of 
view, since its modulus transfer function (MTF) will be 
un-bound and create some aliasing. As we will see, the 
implementation of a diaphragm of appropriate size, mod- 
ifies the slit function in a way which turns out to be ben- 
eficial both from a cross-talk and from a sampling point 
of view. 

Figure [8] shows the BIGRE spaxel conceptually, in- 
dicating its dimensions and the geometrical ray paths. 
The two lenslet-array are implemented as the two sur- 
faces of a single component and the micro-pupil array 
occurs within the component. In principle, it would be 
possible to implement a mask in this micro-pupil image, 
but this option has not been retained in view of com- 
plexity of manufacturing and aligning a system of three 
micro-optical elements (lens, diaphragm, lens). Instead, 
we consider the second lens and the subsequent colli- 
mation optics to be sufficiently large to not significantly 
modify the field transmission, implementing the mask in 
the metapupil image formed onto the spectrograph's dis- 
persion element, see Figure IT2"1 




Input beam. F in 



Fig. 8. — Scheme of a BIGRE spaxel working at the diffraction 
limit with an un-resolved entrance pupil. The first lens lies on a 
focal plane and re-images a micropupil at its focal distance (/i). 
The electric field imaged onto this optical plane is a sine function 
(one dimension) or a Jinc function (two dimensions). This signal 
is filtered by a top-hat transmission function and finally re-imaged 
onto an image plane by the second lens. The distance between 
this intermediate pupil plane and the second lens is its focal length 
(/2). The electric field imaged by this second lens is an un-bound 
signal with upper envelope much steeper than the one proper to a 
sine profile (oc u -1 ) or a Jinc profile (oc u~ 



While the geometrical micropupil size is given by the 
focal ratio of the input beam according to equation (f39|) . 
the characteristic size of the diffractive micropupil is: 



D 



MP 



h 



D s 



(46) 
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In the following, we use a pupil co-ordinate unit, s, 
which is normalized to Dmp, allowing us to discuss the 
size of the pupil diaphragm without worrying about the 
optical design characteristics of the intervening optics. 
For the above assumption concerning relatively undis- 
turbed propagation of the electric field from the microp- 
upils to the spectrograph pupil, we need to ensure that 
the diffracted beam does not get truncated by the second 
microlens edge. For this, a criterion would be to make 
sure the diffractive micropupil is much smaller than the 
spaxel diameter: Dmp << D spaxe i. Plugging this con- 
dition into equation (j46|) . we get the following condition 
on the focal length of the first surface: 



fx « 



D 2 

spaxel 

A ' 



(47) 



(48) 



Introducing a pupil mask defined by: 

PM(s) = U(s/S PM ), 

where II(a;) is a top-hat function with unitary transmis- 
sion within the diameter Spm, which is turn is the pupil 
mask size in units of s. We can express the electric field 
distribution in the exit slit plane as: 



Esiit{u) = FT [e pupt i( s ) ■ PM(s)] 



(49) 



Hence, evoking the convolution theorem and remem- 
bering that the field in the pupil plane is the Fourier 
transform of the field in the spaxel, this can be re- written 
as: 

E slit (u) = n(u/D spaxel ) ® FT [U(s/S P m)} , (50) 

i.e. the convolution between a top-hat function corre- 
sponding to the original spaxel transmission function and 
a Jinc function corresponding to the micropupil mask. 

Finally, the slit function is the square modulus of this 
signal: 

SF(u) = \E s i it (u)\ 2 , (51) 
and the its spectral modulation transfer function is: 



MFT(s) = \FT[SF(u)]\ 



(52) 



8.1. Sampling analysis applied to the BIGRE concept 

According to condition (|19p . the input focal ratio of 
the light coming to the single BIGRE spaxel should be: 



2D 



spaxel 



A, 



(53) 



From the paraxial perspective, input and output fo- 
cal ratios of a BIGRE spaxel are related through the K 
factor as follows: 



F u - 



K 



and the geometric micropupil size returns: 
Dmpg 



A 

F- 



(54) 



(55) 




Fig. 9. — Normalized BIGRE entrance slit function: the example 
shows the case of the IFU optimized for SPHERE in the working 
wavelengths range: 0.95 — 1.35 /im; colors indicate seven distinct 
wavelengths. 



From the diffractive perspective, the output focal ratio 
is fixed only when the size of the pupil mask is fixed 
on the micropupil plane, due to the un-bound nature 
of this micropupil profile. The characteristic size of the 
diffractive micropupil {Dmp) can be parameterized in 
terms of the focal ratio of the first BIGRE lens (F\ ) and 
the spectrograph central wavelength (A c ) as follows: 



D 



MP 



= Spm • A c • Fx . 



(56) 



The diffractive output focal ratio (F ou t) results then 
from the following equation: 



Fout — F~ 



or, exploiting equations ([541 



Dmpg 
Dmp 
1551) and EU) 



D 



spaxel 



K • Spm • A 



(57) 



(58) 



Finally, by equations (|2"Tj) and ([55]) . the spatial cut-off 
frequency of the spectrograph's exit slit becomes: 



K ■ Spm 



'spaxel 



(59) 



Equation (|59| indicates that the actual profile of the 
spectrograph's slit function is no longer a bound signal, 
just because the pupil mask gets a finite size. The actual 
size of the final exit slit function will be then an un-bound 
signal with spatial cut-off frequency depending both on 
the size of this pupil mask, on the de-magnification factor 
of the BIGRE spaxel and the magnification of the re- 
imaging optics. By this analysis, Super-sampling applies 
to the final exit slit function through condition (|2"2")l as 
follows: 



2-d 



pixel 



< 



D 



.-/»/ eel ■ m-IFSj 



K ■ Spm 



(60) 



We can now study the effect of varying the pupil mask 
size on the slit function in term of cross-talk performance 
and on the MTF in term of aliasing. 

Choosing a very large pupil mask, Spm >> 1, cor- 
responds to transmitting the spaxel transmission profile 
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without modification: 
talk is zero. The MTF is a Jinc function with first zero 
at 1.22/s, so sampling this slit with two pixels across 
its width causes aliasing of up to around 15%. On the 
other hand, choosing a very small pupil mask, Spm << 1 
creates a wide slit function with a shape approximately 
equal to an Airy pattern of FWHM ~ D spaxe ij 'Spm- 
The cross-talk is the same as that found for the TIGER 
case, and the MTF is equal to the classical MTF function 
for diffraction limited optical systems. Sampling corre- 
sponding to half of the FWHM is exempt of any aliasing. 

It is somewhat surprising to find in between these two 
extremes, the evolution of the FWHM is not monotonic, 
but passes through a minimum, located at Spm = 2.5. 
At this position, the slit function has a Gaussian-like bell 
shape, and its FWHM is ~ 0.56 • D spaxe i, see Figure fTUl 
The slit function falls off rapidly, and its first secondary 
maximum peaks at values < 0.001. Compared with the 
ones of an Airy function (> 0.01) this ensures a low level 
of cross-talk. The spectral MTF also resembles a Gaus- 
sian function, with a monotonic fall-off, see Figure 111! 
For a sampling of two pixels across the FWHM the alias- 
ing is well below 6%, see Figure [TT] 

The presence of a minimum indicates that the 
size of the slit function could be stable with re- 
spect to variations in wavelength, indicating that the 
pupil mask works chromatically as a pupil apodiza- 
tion (|Jacauinot fc Roisin-Dossierl [1964). This is in- 
deed the indicated in Figure [9j where the slit 
function is plotted for several wavelengths in the range 
0.95 — 1.35 fim. We study the wavelength evolution of 
coherent and incoherent cross-talks in § [9] 




ADIMENSIONAL SIZE OF THE MICROPUPIL MASK (S pu ) 

Fig. 10.— FWHM of the BIGRE entrance slit function profile 
as a function of the pupil mask size Spm- This FWHM gets its 
absolute minimum when Spm = 2.5. Dot-dashed horizontal line 
indicates the asymptotic trend of the FWHM, corresponding to 
the pupil mask sizes towards the limit: Spm — °°- Dashed hor- 
izontal line indicates FWHM value corresponding to the absolute 
minimum Spm = 2.5, this one traced with a dashed vertical line. 



This analysis suggests then that the spectrograph's en- 
trance slit shape can be fixed selecting properly the pupil 
mask's minimum size. Once projected on the final detec- 
tor plane, Super-Sampling can be fixed by imposing that 
two pixels cover the spectrograph's exit slit FWHM: 




i, :ri '-LizEr t-ti-l ft e : ef;~ 



Fig. 11.— MTF of the BIGRE entrance slit: filtering out the 
spatial frequencies above the one corresponding to the entrance slit 
FWHM (dashed horizontal line) produces a limited aliasing error: 
MTF(s) | s > spwhm < 6%. 

where Spwhm is the spatial frequency corresponding to 
a spatial period equal to the slit function FWHM. While, 
according condition ([26]), Hyper-sampling depends on 
the post-coronagraphic pupil size, the spectrograph's 
working wavelengths range and its resolving power. 



TABLE 2 

Independent parameters of a BIGRE-oriented IFS 



Post-coronagraphic pupil size = D 

IFS cut-on wavelength = X m i n 

IFS central wavelength = A c 

Size of the single BIGRE microlens = D spaxe i 

Focal length of the first BIGRE optical surface = fx 

Focal length of the second BIGRE optical surface = fi 

Size of the pupil mask in unit of A c ■ Fx = Smpm 

IFS detector pixel size = d v i xf ,i 

IFS optical magnification = mips 

IFS disperser (2-pixel) resolving power = R 



Finally, as shown in Figure [12l the aim of the optics 
downstream the BIGRE lenslet-array is to re-image the 
entrance slit into the spectrograph's image plane with the 
highest stability and optical quality; for this reason the 
optical design can be fully dioptric. The requested stabil- 
ity is assured imposing the telecentricity of the entrance 
pupil. In turn, this implies that the metapupil forming 
between collimator and re-imaging optics, which is the 
result of the over-position of individual micropupils form- 
ing inside the lenslet-array, has a size equal to the size of 
a single micropupil, once properly magnified by the ratio 
between the equivalent focal length of the collimator op- 
tics and the focal length of the second optical surface of 
the single BIGRE lens. While, the spatial filtering of the 
micropupils is obtained by adopting a unique pupil stop 
placed onto this spectrograph's metapupil plane with a 
physical size (Dps) obtained as follows: 



D 



PS 



D 



MP 



fcoll 



(62) 



2-d 



pixel 



< 



SPWHM 



(61) 



where Dmp is fixed by equation (|56[) . Thus, a suitable 
dispersing device can be inserted in the optical train after 
this pupil stop allowing to image the exit slits as true 
spectra on the spectrograph's image plane. 
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Fig. 12. — BIGRE spectrograph concept: the entrance slits 
plane is filled with the micro-images of the the first surfaces of the 
BIGRE spaxels and the spatial filtering is done in the metapupil 
plane forming between collimator, having focal length fcoll and 
the re-imaging optics, having focal length fcam = miFS ■ fcoll- 

9. BIGRE AND TIGER IFU SOLUTIONS FOR SPHERE IFS 

Coherent and incoherent cross-talks establish the 
actual imaging contrast measured onto the detector 
{Cpixei) with respect to the reference value (C spaxe i) 
depending on the spatial sampling of the post- 
coronagraphic speckle field. Their difference depends on 
cross-talk just because the optical signal imaged by a 
fixed spaxel is spread over a number of detector pixels 
larger than the ones corresponding to it by geometrical 
optical propagation only, in a way which is proportional 
to the levels of cross-talks. When the cross-talk coeffi- 
cients are sufficiently small this difference can be approx- 
imated as: 

Cpixei — C S paxel ~ ~Tl • (CCT + ICT), (63) 

where n is the number of adjacent spaxels around the 
fixed one, while CCT and ICT are the cross-talk coeffi- 
cients defined by equations (|31 [) and (|3"6")) . respectively. 

For the IFS channel of SPHERE the requested cross- 
talk coefficients have been determined through a series 
of simulations devoted to measure the contrast capabili- 
ties of this integral field spectrograph. The result is that 
the impact of cross-talk is well reduced when the Super- 
sampling condition is verified. This fact can be explained 
heuristically remembering the meaning of the cross-talk 
errors over a fixed spectrograph's exit slit: to replace its 
monochromatic intensity with the sum of this intensity 
and the average of the intensities proper to the exit slits 
corresponding to its adjacent spaxels (via the coherent 
cross-talk coefficient) together with the average of the 
intensities proper to the exit slits corresponding to its ad- 
jacent spectra (via the incoherent cross-talk coefficient). 
In the case of Super-sampling, adjacent exit slits do not 
suffer from a mutual shape variations, instead they suffer 
only from mutual differences in intensity due to the input 
post-coronagraphic speckle field. In this way, the resid- 
ual between a fixed exit slit's intensity an its ideal value 
(free from cross-talk errors) becomes small beyond a fixed 
threshold depending on the speckle rejection capabilities 
of the coronagraph. No gain in contrast is then possible 
for further decrements of the cross-talk coefficients. In 
the case of the IFS simulations this threshold returns to 
be 0.01, see Figure [T3l As a conclusion, the IFU solu- 
tions for the IFS of SPHERE should be compliant with 
this specification. 

Figures [U] and [15] show the levels of incoherent cross- 
talks, respectively in the TIGER and BIGRE designs op- 
timized for SPHERE, plotted against wavelength. While 
the incoherent cross-talk is below the 1% threshold for 
both designs, the BIGRE design is clearly superior, show- 
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Fig. 13. — Average 5-cr contrast over an azimuthal area com- 
prised between 0.25 and 1 arc-second from a J = 3.75 (mag) star. 
The simulations are for 10 4 see exposure time and 90 degree field 
rotation. Filled squares are results of the IFS end-to-end simula- 
tions while the solid li ne r epresents the contrast curve expected 
by exploiting equation II63I I in order to obtain C v i xe i via the S- 
SDI calibration technique The adopted post-co ronagraph contrast 
profile is the one presented by [Boccalctti ct al. (2008). 

ing a minimum towards the middle of the range corre- 
sponding to the wavelength at which the pupil mask is 
optimal. The coherent cross-talk is greater than the in- 
coherent one, as expected, but again the BIGRE design 
shows superior performance, and remains well below the 
1% threshold across the spectral range of interest. The 
TIGER design, on the other hand, is not within the spec- 
ified limit. 




0.9 1.0 1.1 1.2 1.3 1.4 

FS WAVELENGTH RANGE = 0.95-1.35 [micron] 



Fig. 14. — Incoherent cross-talk coefficient as a function of the 
wavelength in the range 0.95 — 1.35 fim. Solid line represents the 
BIGRE solution, dashed line the TIGER one and dot-dashed line 
the SPHERE IFS specification. 



Table [3] resumes the solution we found for the BIGRE- 
oriented IFU of SPHERE allowing to reach the requested 
coherent and incoherent cross-talk levels. With this solu- 
tion Hyper-sampling is verified within the whole scientific 
field of view: the Nyquist radius is larger than the radial 
field of view imaged by the spectrograph's optics (1.25 
arc-seconds) . Table ffl] summarizes the solution we found 
for the TIGER-oriented IFU of SPHERE. This one al- 
lows to reach the requested incoherent cross-talk limit 
but not the requested coherent cross-talk limit, while 
Hyper-sampling is well verified as in the previous case. 
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Fig. 15. — Coherent cross-talk coefficient as a function of the 
wavelength in the range 0.95 — 1.35 fim. Solid line represents the 
BIGRE solution, dashed line the TIGER one and dot-dashed line 
the SPHERE IFS specification. 

Based on these results, a BIGRE design is chosen for 
the IFS channel of SPHERE, configured with circular 
spaxels in a hexagonal lattice configuration. 

TABLE 3 

Basic parameters proper to the BIGRE solution op 
SPHERE 



D S paxei = 161.5 fim \ min = 0.95 fim K = 4.1 F out = 12 
miFS = 1-69 d pixel = 18 fim R = 54 



TABLE 4 

Basic parameters proper to the TIGER solution op 
SPHERE 



Dspaxel = 150 fim \ min = 0.95 /im F out = 7 
m IFS = 2.4 d vixel = 18 fim R = 24 




Fig. 16. — Sketch of the selected hexagonal configuration for the 
final spectra of SPHERE IFS. The left-oriented axis is the reference 
on the IFU plane (filled with hexagons representing a portion of the 
spaxels lattice), while the right-oriented axis traces the dispersion 
direction and the black rectangles the spectra imaged onto the 
detector plane. The spectra are 35 pixel long in the dispersion 
direction and the separation to the nearest neighborhood is 5 pixel, 
both in the spectral and in the spatial direction. 



10. COMPARING DIFFERENT BIGRE AND TIGER SPAXEL 
SHAPES AND IFU LATTICE CONFIGURATIONS 

In this Section we compare the slit functions gener- 
ated through the TIGER and BIGRE image propagation, 
computed for different spaxel shapes and lattice config- 
urations of the entire IFU. This comparison is made as- 
suming common spaxel size and wavelength. This analy- 
sis allows to derive the best lenslet-array optical concept 
and the optimum IFU lattice configuration in the ideal 
diffraction limited case, i.e. when the object plane of the 
lenslet-array is an un-resolved entrance pupil. 

The diagnostic quantities exploited for this analysis are 
the amount of coherent and incoherent intensities both 
measured onto the entrance slits plane of the spectro- 
graph, before any chromatical dispersion and re-imaging 
onto a suited detector plane. To this aim, it is important 
to stress the meaning of coherent and incoherent signals 
and the one of their related cross-talk terms. Coherent 
signal is the intensity term due to interference between 
adjacent spaxels measured at any point of the entrance 
slits plane. Such a signal depends on the optical path 
difference between adjacent spaxels only; in this sense 
spaxels can be compared to apertures of a standard grat- 
ing. The coherent cross-talk coefficient is the maximum 
amount of coherent signal, see § 16.11 Incoherent signal is 
the stray intensity terms due to the image propagation 
diffraction effects measured at any point of the entrance 
slits plane. Such a signal depends on the distance be- 
tween adjacent spectra projected onto this plane; in this 
sense this signal depends on the final configuration of the 
spectra on to the detector plane. The incoherent cross- 
talk coefficient is the maximum amount of incoherent 
signal, see ? 16.21 

The comparison between TIGER and BIGRE is per- 
formed for two distinct shapes of the single spaxel (circu- 
lar and square) and for two distinct IFU lattice configu- 
rations (hexagonal and square) . The combination of such 
different shapes and configurations allows to compare the 
TIGER and the BIGRE concepts in term of coherent 
and incoherent signals for standard lenslet-array optical 
setups. It is important to notice that these simulations 
consider as input a normalized signal without amplitude 
and phase differences between adjacent spaxels. In this 
way the results obtained are independent with respect to 
the actual speckle pattern beating the IFU. 

As Figure [T71 indicates . adopting a bi-logarithmic scale, 
the single BIGRE slit gets an intensity profile steeper 
than the one proper to the single TIGER slit both in 
the case of circular and square shapes. More in detail, 
the upper envelope to the slit intensity profile proper to 
a circular TIGER-oriented spaxel is a power law with 
index equal to —3, while the same quantity for a square 
TIGER-oriented spaxel is a power law with index equal 
to —2 along the aperture side and with index equal to —4 
along its diagonal. At contrary, the upper envelope of the 
slit intensity profile proper to a circular BIGRE-oriented 
spaxel is not a power law (only its asymptotic tail is fitted 
quite well with a power law having index oc —4.5); the 
same quantity is not a power law in the case of a square 
BIGRE-oriented spaxel too (only its asymptotic tail is 
fitted quite well with a power law with index ex —3 in 
the direction of the aperture side and index cx —6 along 
its diagonal). 
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Fig. 17. — Normalized BIGRE and TIGER, slit functions comparison. The adopted lenslet aperture is always 161.5 /im i.e. the one 
proper to the BIGRE solution for SPHERE (see Tabled; the adopted wavelength is always 1 (im. First column shows the images given by 
individual spaxels in a bi-logarithmic scale. Second column shows the relative profiles (in the case of square spaxel shape diagonal profiles 
are presented in green color) . The power laws fitting the slit functions proper to a cicular and a square TIGER-oriented spaxel are indicates 
for reference by solid lines. These simulations consider Fraunhofer propagation only and assume as input signal of the lenslet-array an 
un-resolved entrance pupil. Finally, results are independent to the detector pixel scale just because images refer to the spectrograph's 
entrance slits plane. 
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CIRCULAR BIGRE SPAXEL WITHIN A HEXAGONAL IFU LATTIC 
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CIRCULAR TIGER SEAXEL WITHIN A HEXAGONAL IFU LATTIC 
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Fig. 18. — Normalized BIGRE and TIGER, slit functions comparison. The adopted lenslet aperture is always 161.5 /im i.e. the one 
proper to the BIGRE solution for SPHERE (see Tabled; the adopted wavelength is always 1 (im. First column shows the images given by 
an individual spaxel together with its adjacent ones and for different IFU lattice configurations. Second column shows the relative profiles: 
black color indicates the slit function (the signal represents the detection expected if only the central spaxel were illuminated); red color 
indicates the signal due to adjacent spaxels only; blue color is the coherent signal arising from the interference of the signal proper to 
the central spaxel with all the adjacent ones. These simulations consider Fraunhofer propagation only and assume as input signal of the 
lenslet-array an un-resolved entrance pupil. Finally, results are independent to the detector pixel scale just because images refer to the 
spectrograph's entrance slits plane. 
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The result is that the BIGRE-oriented circular aper- 
ture within a hexagonal lattice configuration allows a 
superior suppression of coherent and incoherent signals, 
while the slits generated by a circular TIGER-oriented 
aperture in a hexagonal lattice are similar — in this 
context — to the ones generated by a square BIGRE- 
oriented aperture in a square lattice. Finally, the slits 
generated by a square TIGER-oriented aperture in a 
square lattice are the worst in term of coherent and inco- 
herent signals suppression, see Figm-efTS] Hence, the con- 
tribution of non-adjacent spaxels can be neglected when 
evaluating the cross-talk signals in the case of a BIGRE 
spectrograph, just because the power laws fitting — in 
a bi-logarithmic plot — the intensity distribution proper 
to the TIGER slit functions do not fit at all the one 
proper to the BIGRE slit functions. At contrary, the in- 
tensity distribution proper to the BIGRE slit functions 
can be only approximated with lower index power laws. 
Thus, what for a TIGER lenslet-array represents an es- 
timate only, for a BIGRE lenslet-array it gives realistic 
measures of the signals due to the spectrograph's slit 
functions cross-talk. 

11. CONCLUSIONS 

By integral field spectroscopy it is possible to realize 
the S-SDI ca l ibrati on technique in the way proposed by 
iBerton et al.l (|2006[ ). and — at least in a few cases — to 
get the spectrum of candidate extrasolar giant planets 
adopting suited s pectral de- co n voluti on recipes, as the 
one proposed by iThatte et al.l (|2007f ). However, these 
techniques can increase the contrast performances only 
when several sampling conditions, both in the spatial and 
in the spectral domain of the speckle field, are verified. 

In this context, our effort has been to discuss in gen- 
eral terms the critical sampling conditions needed to deal 
with a speckle field data cube before applying on it the S- 
SDI calibration technique or any spectral de-convolution 
recipe. To this purpose, we evaluated the impact of the 



cross-talk as function of various parameters of a lenslet- 
based integral field spectrograph, especially in the case 
of trying to minimize the number detector pixels (which 
is an issue in general for IFS) in the case of strong specifi- 
cations, as the ones requested for high-contrast imaging. 
For this reason we conceived a new optical scheme — we 
named BIGRE — and characterized it in the specific case 
of the IFS channel foreseen inside SPHERE, showing that 
a BIGRE-oriented spectrograph is conceptually feasible 
by standard dioptric optical devices. Once applied to 
the technical specifications of this instrument, a BIGRE 
integral field unit is able to take into account the effects 
appearing if a lenslet-array is used in diffraction-limited 
conditions. Specifically, we proved here that coherent 
and incoherent cross-talk coefficients reach values deeper 
than for a TIGER IFU when applied to the same optical 
frame. More in general, the comparison between the BI- 
GRE and the TIGER spaxel concept has been pursued in 
terms of coherent and incoherent cross-talk suppression, 
adopting a common size for the single aperture and a 
fixed monochromatic wavelength for the wavefront prop- 
agation. In the ideal case of uniform illumination with 
un-resolved entrance pupil, the circular BIGRE spaxel 
within an hexagonal IFU lattice configuration shows to 
be the optimal solution among the ones we investigated. 
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